APS/123-QED 



Magnons and electromagnons in a spin-lattice-coupled frustrated magnet CuFe02 

as seen via inelastic neutron scattering 

Taro NakajimajB Azusa Sunolll and Setsuo Mitsuda 
Department of Physics, Faculty of Science, Tokyo University of Science, Tokyo 162-8601, Japan 

Noriki Terada 

National Institute for Materials Science, Tsukuba, Ibaraki 305-0044, Japan 



Shojiro Kimura 

Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan. 

Koji Kaneko and Hiroki Yamauchi 

Quantum Beam Science Directorate, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan 

We have investigated spin-wave excitations in a four-sublattice (4SL) magnetic ground state of 
a frustrated magnet CuFe02, in which ' elect romagnon' (electric- field-active magnon) excitation 
has been discovered by recent terahertz time-domain spectroscopy [Seki et al. Phys. Rev. Lett. 
105 097207 (2010)]. In previous study, we have identified two spin- wave branches in the 4SL 
phase by means of inelastic neutron scattering measurements under applied uniaxial pressure. [T. 
Nakajima et al. J. Phys. Soc. Jpn. 80 014714 (2011) ] In the present study, we have performed 
high-energy-resolution inelastic neutron scattering measurements in the 4SL phase, resolving fine 
structures of the lower-energy spin- wave branch near the zone center. Taking account of the spin- 
driven lattice distortions in the 4SL phase, we have developed a model Hamiltonian to describe 
the spin- wave excitations. The determined Hamiltonian parameters have successfully reproduced 
the spin-wave dispersion relations and intensity maps obtained in the inelastic neutron scattering 
measurements. The results of the spin- wave analysis have also revealed physical pictures of the 
magnon and electromagnon modes in the 4SL phase, suggesting that collinear and noncollinear 
characters of the two spin- wave modes are the keys to understand the dynamical coupling between 
the spins and electric dipole moments in this system. 

PACS numbers: 75.30.Ds, 78.70.Nx, 75.80.+q, 75.85.+t 



I. INTRODUCTION 

A triangular lattice antiferromagnet CuFe02 (CFG) 
has recently attracted increasing attention due to its 
cross-correlated phenomena arising from spin frustration. 
CFG has a delafossite structure, in which triangular lat- 
tice layers of magnetic Fe^+ ions are separated by non- 
magnetic G^~-Cu+-G^~ dumbbells. From the electronic 
state of Fe^+ ions {S = 5/2, L = 0) and antiferromag- 
netic interactions between them, one might expect a non- 
collinear three-sublattice magnetic ground state so called 
"120°-structure".^ However, CFG exhibits, in a ground 
state, a collinear four-sublattice (4SL) antiferromagnetic 
order in which the magnetic moments of Fe^+ ions are 
confined along the c axis^ (see Fig. [IJa)). Recent theo- 
retical studies have pointed out that in frustrated mag- 
nets, collinear magnetic orderings can be stablized by 
strong spin-lattice coupling. Actually, synchrotron ra- 
diation x-ray diffraction studies by Terada et al.^ and 
Ye et al.^ have revealed that CFG exhibits 'spin-driven' 
crystal structural transitions; while the crystal structure 
of CFG is a trigonal structure (space group R3m) in the 
paramagnetic (PM) phase, it turns to be a monoclinic 
structure in the 4SL phase. This implies that the lattice 
degree of freedom is the key to stabilize the 4SL magnetic 
order. 



CFG is also known as a spin-driven magneto-electric 
(ME) multiferroic,^ in which a screw-type magnetic order 
induced by nonmagnetic substitution or application of 
magnetic field breaks inversion symmetry of the system, 
and triggers ferroelectricity^^ — Furthermore, recent ter- 
ahertz time-domain spectroscopy has discovered 'electro- 
magnon (electric-field-active magnon)' excitation in the 
nonferroelectric 4SL phase, indicating a dynamical cou- 
pling between spins and electric dipole moments in this 
system. 

To understand these various cross-correlated phenom- 
ena, it is fundamental to determine the magnetic interac- 
tion parameters of this system. Ye et al. have performed 
inelastic neutron scattering measurements to determine 
these parameters from spin-wave dispersion relations in 
the 4SL phase. However, in the previous work, there 
is some ambiguity in the identification of the spin-wave 
branches because of a problem arising from magnetic 
domain structures in this system. Specifically, in the 
4SL phase, CFG has three types of magnetic domains 
whose magnetic propagation wave vectors are described 

as (i,i,f), (-^,i,f) il^-h^) ^si^g a hexago- 
nal basis. These wave vectors are crystallographically 
equivalent to each other because of the threefold rota- 
tional symmetry of the original trigonal crystal structure. 
Therefore, the magnetic excitation spectrum reported in 
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the previous stud}4i^ is a mixture of spin- wave spectrums 
corresponding to the three different orientations of the 
magnetic domains. 

Quite recently, we have demonstrated that a 'single- 
domain' 4SL phase can be realized by uniaxial pressure 
applied perpendicular to the c axis.^^ We performed in- 
elastic neutron scattering measurements under applied 
uniaxial pressure, and identified two spin- wave branches 
in the single-domain 4SL phase. We have also pointed 
out that the lower-energy branch splits into two branches 
near the zone center. However, the wave- vector (Q) de- 
pendence of the splitting was not fully investigated. 

In the present study, we have thus performed high- 
energy-resolution inelastic neutron scattering measure- 
ments to resolve the fine structure of the spin-wave dis- 
persion relations in the 4SL phase. Taking into account 
of the spin-driven lattice distortion, we have developed 
a model Hamiltonian to describe the observed spin-wave 
spectrums. The results of the spin- wave calculation have 
revealed the microscopic picture of the magnon and elec- 
tromagnon modes in the 4SL phase. 



II. EXPERIMENT AND RESULTS 

We used the single-crystal CFO sample identical to 
the sample used in Ref. 17. The sample has been set 
in a uniaxial pressure cell developed by Aso et al}^ An 
uniaxial pressure of 10 MPa has been applied to the [110] 
surfaces of the sample in order to produce the single- 
domain 4SL state in which the magnetic domains with 
the wave vector of (|, |, f ) dominate over the others.— 

We used a cold neutron triple axis spectrometer 
LTAS(C2-1) installed at JRR-3 in Japan Atomic En- 
ergy Agency, Tokai, Japan. The energy of the scat- 
tered neutrons was fixed to Ef = 2.5 meV. A horizontal 
focusing analyzer was employed. Energy resolution at 
elastic position was 0.06 meV (full width at half maxi- 
mum). The higher-order contaminations were removed 
by a cooled Be-filter placed in front of the sample and 
a room-temperature Be-filter placed in front of the ana- 
lyzer. 

As mentioned in the introduction, we have already 
identified two spin- wave branches in our previous work.^^ 
We have referred to the higher and lower energy branches 
as the HE- and LE-branches, respectively. To resolve the 
fine structure of these branches, we performed constant- 
Q scans at several reciprocal lattice points at T = 2 
K in the 4SL phase, as shown in Figs. [Hdj-IHg). We 
found that the scattering profiles apparently have two- 
peak structures around the energies corresponding to the 
LE-branch. On the other hand, splitting of the HE- 
branch was not observed within the present experimental 
resolution, as shown in the inset of Fig. [Hf). Therefore, 
we focused on the LE-branch, and systematically carried 
out constant-Q scans along |) line in the range of 

0.15 <H < 0.35 and < E < 2.0 meV. In Fig. [U^c), we 
have plotted the peak positions of the observed scattering 
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FIG. 1: (Color Online) (a) A magnetic unit cell of the 4SL 
magnetic structure, am, hm and Cm denote the monoclinic ba- 
sis, (b) The reciprocal lattice map of the (if, if, L) scattering 
plane in CFO. (c) The peak positions of the constant-Q scans 
combined with the results in Ref. ITtI . The horizontal bars de- 
note full widths at half- maximums of the scattering profiles 
in the previous study, [(d)- (g)] The scattering profiles of the 
constant-Q scans at (d) (0.18,0.18,1.5), (e) (0.21,0.21,1.5), (f) 
(0.25,0.25,1.5) and (g) (0.29,0.29,1.5), at T = 2 K in the 4SL 
phase. 



profiles together with the results in Ref. flTl, revealing the 
dispersion curves of the two split branches. Hereafter, we 
refer to the split branches having lower and higher ener- 
gies as "LEI" and "LE2" branches, respectively. 



III. CALCULATIONS AND DISCUSSIONS 

A. Model Hamiltonian for the Spin- Wave 
Excitations in the 4SL Phase 

The first calculation for the spin- wave dispersion rela- 
tions in the 4SL phase has been presented by Fishman^i^ 



Trigonal 
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FIG. 2: (Color Online) [(a)- (h)] Paths of the exchange inter- 
actions [(a)- (d)] in the trigonal symmetry and [(e)- (h)] the 
monoclinic symmetry of the crystal structure of CFO. 



Although the calculation, in which three in-plane ex- 
change interactions (Ji, J2 and J3), a exchange interac- 
tion between adjacent layers (J^) and a uniaxial single 
ion anisotropy {D) are employed, successfully reproduces 
the dispersion relation of the LE branch, it does not that 
of the HE branch. Specifically, in Ref. 19, the exci- 
tation bandwidth for the HE branch along (i^,i^, 3/2) 
line is calculated to be ~ 2.5 meV (from 2.5 to 5.0 
meV), while that has been determined to be ~ 1.2 meV 
(from 2.0 to 3.2 meV) in our experiment. In subse- 
quent study, Fishman et al. have corrected the paths 
of the exchange interactions between the adjacent layers, 
and have presented revised calculations for the spin- wave 
dispersion relations. However, the dispersion relation 
of the HE branch has not been discussed. These previ- 
ous calculations^^'^^ are based on the assumption that 
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FIG. 3: (Color Online) (a) Isosceles and (b) scalene triangle 
lattices with the 4SL magnetic order. 



all of the exchange interactions are isotropic, as shown in 
Figs. [2fa)-|2fd). In other words, effects of the spin-driven 
lattice distortions on the exchange interactions are not 
taken into account. 

On the other hand, Kimura et al. have recently 
discussed the effects of lattice distortions on ESR sig- 
nals, i.e., spin- wave energies at zone center, in the 4SL 
phase^^l They have performed multi-frequency ESR mea- 
surements in the 4SL phase, and have analyzed the re- 
sults using the 'scalene triangle model' proposed by the 
synchrotron radiation x-ray diffraction study by Terada 
et al^ In this model, the nearest neighbor exchange inter- 
actions, Ji, splits into three inequivalent interactions due 
to the monoclinic lattice distortion and displacements of 
the 0^~ ions, as shown in Fig. [3l^b). These deformations 
are expected to enhance the antiferromagnetic exchange 
interactions between neighboring up and down spins, and 
on the contrary, to reduce those between neighboring two 
up (or down) spins, as was discussed in Refs. and [22I . 
Consequently, the splitting of Ji contributes to stabilize 
the 4SL magnetic ordering. 

Quite recently, Kimura et al have calculated spin- wave 
dispersion relations in the 4SL phase using the scalene tri- 
angle model, assuming that Ji splits into three inequiv- 
alent interactions, while J2, J3 and Jz are isotropic 
They have shown that the energy bandwidth of the HE 
branch along (i^,i^, 3/2) line is calculated to be ~ 1.0 
meV (from 1.8 to 2.8 meV). In addition, their calcula- 
tions have qualitatively reproduced overall features of the 
results of the inelastic neutron scattering measurements 
in the multi-domain state. However, there still remain 
some discrepancies between the 'shapes' of the calculated 
and observed dispersion curves for the single-domain 4SL 
phase. 

In the present study, we start from three in-plane ex- 
change interactions Ji , J2 and J3 and a inter-layer ex- 
change interaction J^. First, we introduce effects of the 
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JiS 






J2S 


J2S 






-0.455 


-0.422 


-0.150 


-0.100 


-0.106 


-0.338 


-0.372 


JzS 


J'zS 


DS 


ES 








-0.187 


-0.167 


0.160 


-0.035 









TABLE I: The Hamiltonian parameters (in meV). 



trigonal to monoclinic structural transition on all of the 
exchange interactions. As shown in Figs. [2fe)-|2fh), the 
monoclinic symmetry of the crystal allows the exchange 
interaction {i — l,2,3,z) to split into Ji and J[. Sec- 
ond, we have applied the scalene triangle model to the 
nearest neighbor exchange interactions, and therefore J( 
splits into J{ and J'{^ as shown in Figs. [3l^a) and[3fb). 
Although the scalene triangle distortion might affect the 
other distant interactions, we have assumed that it af- 
fects only the nearest neighbor interactions in the present 
analysis. 

As for the single ion anisotropy, Fishman et al. have 
employed only a uniaxial anisotropy, 1^.16,19,20 q^^ ^^le 
other hand, Kimura et al. have argued that the in-plane 
anisotropy E is necessary to explain the results of their 
ESR measurements.^^ We have thus employed D and E 
terms in the present analysis. Therefore, we have written 
a Hamiltonian for the 4SL phase as 



3 3 

i,j n=l i,j n=l 

-Y,D{Stf-J2E[{S^f-{S!)% (1) 

where the x, y and z axes are defined to be parallel to 
[110], [110] and [001] directions of the crystal, respec- 
tively. 

B. Calculation for the Spin- Wave Spectrum 

We have calculated the spin-wave energies in the 
4SL phase, applying a Holstein-Primakoff 1/S expansion 
about the classical limit to the Hamiltonian of Eq. ([1]). 
Along the previous work by Fishman, we refer to the 
four magnetic sublattices as a, /3, 7 and J, as shown in 
Fig. [3l^b). We express the spins Si on the sublattice a 
(/3,7 and S) using the boson operators ai {Pi^ji and di). 

The procedures of the calculation for the spin- wave en- 
ergies are essentially the same as those in Ref. To'. The 
spin- wave energies at wave vector e^, are obtained by 
solving the Heisenberg equations of motion for the vec- 
tor Vq = {aq,Pq,-fl,Sl,alq,plq,-f-q,5-q), whcrC tt^, Pq 
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FIG. 4: (Color Online) Spin-wave dispersion relations along 
the (H,H, |) line calculated with the parameters in Table |I] 
and the observed spin- wave energies. 



and so on are the Fourier-transformed boson operators. 
The equation of motion for Vq can be written using the 
8x8 matrix M{q) as idvq/dt = —[H,Vq] = M{q)vq. 
Diagonalizing the matrix M{q)^ we have obtained the 
spin-wave energies in the 4SL phase. 

To find a suitable set of magnetic interaction parame- 
ters, we have prepared 5 x 10^ sets of the parameters by 
Monte Carlo method, and have calculated spin-wave 
energies e^., where i stands for a observed data point. 
For each set of the parameters, we have calculated a sum 
of squared residuals = ^ii^qi ~ ^gj^? where Eq. is a 
observed spin- wave energy, and then have searched a set 
of the parameters having a small value of - 1^ should 
be noted that in principle, a finite in-plane anisotropy 
[E term) splits not only the LE branch but also the HE 
branch into two branches. Therefore, we have four cal- 
culated spin- wave branches referred to as LE1-, LE2-, 
HEl- and HE2-branches. Since we could not observed 
the splitting of the HE-branch in our inelastic neutron 
scattering measurements, we have used Eq. of the ob- 
served HE-branch for both of the calculated HEl- and 
HE2-branches, in the evaluation of . 

Examining the 5 x 10^ sets of the parameters, we have 
obtained a set of parameters which successfully repro- 
duces the observed dispersion relations. Table U shows 
the values of the Hamiltonian parameters. In Fig. HI we 
show the calculated spin- wave dispersion relations along 

3/2) line with the observed spin-wave energies. 
Using these parameters, we have also calculated inelastic 
neutron scattering cross-sections described as 
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FIG. 5: (Color Online) Observed intensity maps of the inelastic neutron scattering measurements along (a) the {H, if, 3/2) and 
(c) (1/4, 1/4:, L) lines (taken from Ref. JJi). Calculated intensity maps along (b) the (if, if, 3/2) and (d) (1/4, 1/4, L) lines for 
a single domain 4SL phase, (e) Spin-wave dispersion relations for a magnetically multi-domain 4SL phase calculated from the 
parameters in Tabled (f) Calculated intensity map for the multi-domain state in which volume fractions of the three domains 
are assumed to be equal to each other, (g) The reciprocal lattice map of the (if, if, L) scattering plane. Arrows denote the 
directions of the (if, if, 3/2) and (1/4, 1/4, L) lines. 



where, k is the scattering vector defined SiS k, = kf — ki. 
kf and ki are wave vectors of the scattered and incident 
neutrons, respectively, /(a^) is the magnetic form factor 
of a Fe^+ ion.^^ f<iz is the z-axis component of an unit vec- 
tor of the scattering vector. R denotes a position of a spin 
of a Fe^+ ions. To obtain resolution-convoluted inelas- 
tic neutron scattering spectrums, we have employed the 
Cooper-Nathans type resolution function*^ As shown in 
Figs. [3a)-|5l(d), we have demonstrated that the calcu- 
lated intensity maps show good agreement with the ob- 
served data (taken from Ref. IITI ). 

It should be noted that the spin-wave dispersion rela- 
tions in the (H^H^L) plane does not depend on J2, be- 
cause the direction of J2 is perpendicular to the (if, if, L) 
plane. To determine J2, we have calculated the spin- wave 
spectrums in the magnetically multi-domain state. Be- 
cause the directions of J2 in magnetic domains with the 
wave vectors of (— ^, |, |) and (|, — ^, f ) are not perpen- 
dicular to the (if, if, L) scattering plane, the spin- wave 
spectrums belonging to these domains contain informa- 
tion on J2. We have adjusted the value of J2 so that 
the calculated multi-domain spectrum is close to the ob- 
served data in Ref. As a result, J2S is determined 
to be —0.106 meV. In Figs [Sfe) and[5ff), we show the 
calculated spin- wave dispersion relations and a intensity 
map for an inelastic neutron scattering measurement in 
the multi-domain state in which the volume fractions of 



the three domains are assumed to be equal to each other, 
respectively. 

We now discuss the determined values of the Hamilto- 
nian parameters. Among the nearest neighbor exchange 
interactions, the magnitude of the exchange interactions 
connecting two up (or down) spins, | J{'|, is determined to 
be smaller than | Ji | and | J[ \ which connect up and down 
spins. This is consistent with the nature of the scalene 
triangle distortion; specifically, the nearest neighbor ex- 
change interactions vary so as to lower the total exchange 
energy of the 4SL magnetic structure. It is worth men- 
tioning that previous x-ray diffraction studies on CFO^^- 
and CuFei_a;Ala^02^^ have suggested that this scalene 
triangle distortion occurs not only in the 4SL phase, but 
also in the screw- type helimagnetic (multiferroic) phase. 
Haraldsen et al. have recently presented an inelastic neu- 
tron scattering and theoretical study on magnetic exci- 
tations in the helimagnetic phase of Ga-doped CFO.^^ 
They have assumed that the oxygen displacements also 
break the equilateral symmetry of Ji so as to lower the 
total exchange energy, in the helimagnetic phase. Their 
calculations have successfully reproduced magnetic exci- 
tations from the noncollinear magnetic structure. This 
ensures that the oxygen displacements significantly affect 
the competing magnetic interactions in this system. 

We have also found that the second and third neigh- 
bor in-plane exchange interactions and inter-plane inter- 
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FIG. 6: (Color Online) Schematic drawing of the motion of 
the spins in (a) the LE- and (b) HE-modes. (c) Calculated 
electric polarization in a magnetic unit cell for the LE- and 
HE-modes. Inset shows positions of the four inequivalent oxy- 
gen sites in a magnetic unit cell. [110]-plane-projections of the 
spins of the (d) LE- and (e) HE-modes. 



actions have finite anisotropy, specifically, J2/J2 = 1-05, 
jyjs = 1.10 and J'JJ^ = 0.89. This suggests that the 
effects of the monoclinic structural transition, which are 
not taken into account in the previous studies, ^^li^Ti^ are 
important to explain the spin- wave excitation in the 4SL 
phase. 



C. Microscopic Picture of the "Electromagnon" 
Excitation 

As mentioned in introduction, Seki et al have re- 
cently discovered the "electromagnon" excitation in the 
4SL phase. Specifically, their optical spectroscopy have 
detected two groups of signals corresponding to the zone- 
center spin-wave modes belonging to the LE- and HE- 
branches. Hereafter, we refer to these spin- wave modes as 
LE- and HE-modes, respectively. Measuring these signals 
with various light polarization configurations, they have 
concluded that the HE- mode is "electromagnon" driven 
by ac electric field, and the LE-mode is conventional 
" magnon" . To investigate the mechanism of the " electro- 
magnon" excitation, it is indispensable to establish phys- 
ical pictures of these spin- wave modes. In the following 
discussion, we have neglected the in-plane anisotropy, 
for simplicity, and therefore there are only two spin- wave 
modes corresponding to the LE- and HE-modes, because 
the LEI- and LE2-modes (the HEl- and HE2-modes) de- 
generate without the E term. 

In Figs. [6fa) and[6l^b), we have shown the schematic 
drawings of the LE- and HE-modes at zone center, which 
are obtained by the spin-wave analysis without the E 
term. These pictures are consistent with the spin-wave 



modes presented in the recent ESR study by Kimura 
et al?^ Each spin rotates about the z axis with angu- 

LE/HE LE/HE 

lar frequency of cJq , where cJq is the spin- wave 
frequency of the LE/HE-branch at (7 = 0. In the LE- 
mode, the up spins (or the down spins) rotate keeping a 
"collinear" configuration. On the other hand, in the HE- 
mode, the two up spins (or two down spins) are coupled 
in a " noncollinear" configuration, in a magnetic unit cell. 
In Figs. El^d) and[6fe), we show [110]-plane-projections 
of the spins depicted in Figs. [6l^a) and[6jb), respectively. 
We found that the HE-mode has a screw-like spin tex- 
ture but the LE-mode does not, as was predicted in Ref. 
[15I . Moreover, the spin-helicity of the screw-like spin tex- 
ture is reversed by rotating the spins about the z axis by 
TT. Keeping the relationship between the ferroelectricity 
and the screw-type magnetic ordering in the multifer- 
roic phase of CFO^^^« in mind, we anticipate that the 
noncollinear nature of the HE-mode accounts for the dy- 
namical ME-coupling in the 4SL phase. 

To quantitatively investigate the difference between 
the LE- and HE-modes, we have applied, to the two 
modes, the d-p hybridization model, which successfully 
explains the ferroelectricity in the screw-type helimag- 
netic phase of this system. This model requires three 
noncollinear spins surrounding an oxygen ion as a mini- 
mal unit to produce uniform electric polarization. Hence, 
we have calculated a sum of electric polarizations at four 
inequivalent oxygen sites, each of which is surrounded 
by three Fe ions, in a magnetic unit cell (see inset of 
Figini^c)). The procedures of the calculation are essen- 
tially the same as those in Refs. [30| and [32|. Fig. [6fc) 
shows the calculated electric polarizations along the [110] 

LE /HE 

direction as functions of cJq t for the LE- and HE- 
mode, revealing that the HE-mode generates electric po- 
larization oscillating with the same angular frequency as 
cJq^^, but the LE-mode does not. This oscillating electric 
polarization can be coupled with the ac electric fields. 
We have thus concluded that this must be the micro- 
scopic mechanism of the electromagnon excitation in this 
system. 

IV. CONCLUSION 

We have investigated spin-wave excitations in the 4SL 
magnetic ground state of a frustrated magnet CuFe02. 
By high-energy-resolution neutron inelastic scattering 
measurements, we have revealed the fine structure of the 
spin- wave dispersion relations in the 4SL phase. Taking 
account of the trigonal to monoclinic structural transi- 
tion and the scalene triangle distortion in the 4SL phase, 
we have developed a model Hamiltonian, and have deter- 
mined the Hamiltonian parameters. We have successfully 
reproduced the spin- wave dispersion relations and the in- 
tensity maps for the inelastic neutron scattering measure- 
ments in the single-domain and multi-domain 4SL phase. 
Consequently, we have concluded that the the spin-driven 
lattice distortions significantly affect the exchange inter- 
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actions in this system. 

Using the results of the spin- wave analysis, we have 
presented the physical pictures of the zone-center spin- 
wave modes belonging to the LE- and HE-branches, and 
have discussed the mechanism of the recently discovered 
electromagnon excitation. The pictures of the LE- and 
HE-modes suggest that the symmetry of the spin-wave 
modes are the key to understand this kind of dynamical 
ME coupling. Applying the d-p hybridization model^^ 
to the two modes, we have demonstrated that only the 
HE-mode generates the oscillating electric polarization, 
which can be coupled with ac electric field, but the LE- 
mode does not. This indicates that the d-p hybridization 
model, in which spins and electric dipole moments are 



mediated by a spin-orbit coupling, can explain not only 
the static coupling between the magnetic ordering and 
the ferroelectricity, but also dynamical ME-coupling, i.e., 
electromagnon excitation. 
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